Shear-driven size segregation of granular materials: modeling and experiment 
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Granular materials segregate by size under shear, and the ability to quantitatively predict the 
time required to achieve complete segregation is a key test of our understanding of the segregation 
process. In this paper, we apply the Gray-Thornton model of segregation (developed for linear 
shear profiles) to a granular flow with an exponential profile, and evaluate its ability to describe 
the observed segregation dynamics. Our experiment is conducted in an annular Couette cell with a 
moving lower boundary. The granular material is initially prepared in an unstable configuration with 
a layer of small particles above a layer of large particles. Under shear, the sample mixes and then 
re-segregates so that the large particles are located in the top half of the system in the final state. 
During this segregation process, we measure the velocity profile and use the resulting exponential 
fit as input parameters to the model. To make a direct comparison between the continuum model 
and the observed segregation dynamics, we locally map the measured height of the experimental 
sample (which indicates the degree of segregation) to the local packing density. We observe that the 
model successfully captures the presence of a fast mixing process and relatively slower re-segregation 
process, but the model predicts a finite re-segregation time, while in the experiment re-segregation 
occurs only exponentially in time. 

PACS numbers: 45.70.Mg, 47.57.Gc, 81.05.Rm, 64.60.ah 



I. INTRODUCTION 

Granular materials have long been known to segre- 
gate by size, shape, density, and other material proper- 
ties, whether driven by shear or vibration [1, 2, 3J. The 
case of size-segregation in shear flow is particularly im- 
portant, as it arises in such diverse situations as indus- 
trial chute flows, rock avalanches, and rotating tumblers 
[4j |5j |6j |7j |8] . Under such shear, large particles typically 
rise and small particles descend. The dominant mech- 
anism in such cases is thought to be percolation-based, 
where the granular flow acts as a sieve through which the 
small particles preferentially fall, but arguments based on 
kinetic theory have been proposed as well [13 [TU] . A num- 
ber of groups |TTJ [TU [T3] have developed continuum mod- 
els to describe vertical size-segregation within the con- 
stant shear-rate flows typical of free-surface avalanches. 
The aim of this paper is to apply a recent modification 
[H] of the Gray-Thornton model [Jj|J to a boundary- 
driven granular system where the shear rate is a nonlin- 
ear function of depth. We quantitatively investigate the 
model's ability to capture experimentally-observed seg- 
regation dynamics and the associated changes in overall 
packing density. 

We perform experiments on binary granular mixtures 
of spherical particles confined in an annular Couette ge- 
ometry (see Fig. [IJ. The system starts in an unstably- 
stratified state with small particles over large particles 
and progresses to the reverse, stable, configuration as the 
lower boundary is rotated. Key advantages of this geom- 
etry are the ability to run continuously without the need 
to feed material, as would be the case for chute flows, 
and the ability to start and end the experiment in well- 
defined states. As is commonly the case for sheared gran- 



ular materials [TS] , we observe that the velocity profile is 
not a linear function of depth as assumed in [T2l[T3l[Tfj] . 
but instead decays exponentially away from the shearing 
surface. Correspondingly, the segregation rate is not uni- 
form, but is higher near the bottom rotating plate. We 
use high-speed digital imaging and particle-tracking to 
measure the velocity profile u(z) in the experiment and 
use the resulting fit to an exponential form as input to 
the model. 

Kinetic sieving models of segregation [TTJ [TSJ [T3] are 
based on the notion that vertical size-segregation is prin- 
cipally caused by small particles preferentially falling into 
gaps created by the relative motion (shear) of the par- 
ticles beneath them. Secondarily, the large particles are 
forced upward in a process sometimes called squeeze ex- 
pulsion. In particular, Bridgwater et al. [11] argue that 
the vertical velocity of small particles due to segrega- 
tion (the segregation velocity) should be proportional to 
the shear rate 7 = \du/dz\ as well as some function of 
the concentration <p of small particles. In the context of 
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FIG. 1: Schematic of experimental apparatus (not to scale) 
showing initial configuration and coordinate system. 
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chute flow and avalanches, both Savage and Lun [H] and 
Gray and Thornton 1 13] assume a segregation velocity 
which is additionally proportional to the concentration 
of large particles, (1 — </)), since large particles provide 
the gaps into which the small particles fall. The choice 
of a constant of proportionality (known as the segregation 
rate) comes from an assumption of constant shear rate, 
whereas for many important granular systems it instead 
falls off exponentially [TS] . 

Since the local shear rate is set by the local horizontal 
velocity, we measure u(z) directly from the experiment. 
In our measurements of the horizontal velocity profile, 
described in Sjll] we observe that it falls off exponentially 
away from the shearing surface, is the same for both par- 
ticle sizes, and is approximately time-independent except 
for an initial transient. Using the steady-state velocity 
profile u{z) to determine the shear rate, we calculate an 
exact solution, described in May et al. [14] and summa- 
rized in pil| The concentration of small particles, calcu- 
lated as a function of depth and time, predicts a rarefac- 
tion wave that mixes the particles over short time. As 
small particles reach the bottom plate, a growing layer 
of small particles is formed, bounded by a shock wave 
(a discontinuity in concentration) propagating upwards. 
Correspondingly, large particles eventually reach the top 
of the annulus, and a second shock propagates down- 
ward. The two shocks meet in finite time, after which 
the material is completely re-segregated. 

While it is not possible to monitor the local concentra- 
tion field <fi(z,t) within the experiment, we measure the 
progression of the mixing and segregation processes via 
the compaction and subsequent re-expansion of the ag- 
gregate. To relate the dynamics of the exact solution to 
the experiment, we postulate that the local concentration 
(f) of small particles determines the local packing density 
p, as described for static packings in Kristianscn et al. 
|17j . This concentration map p(4>) allows us to model the 
change in the measured height H(t) of the aggregate as a 
function of time. We compare this proxy height H(t) to 
the time evolution of the experimentally measured height 

in drv 



In jjv] we evaluate the successes and failures of the 
model. To make a quantitative comparison, we set two 
of the three free model parameters from the observed 
system height and the magnitude of the transient com- 
paction. The third parameter is the constant of propor- 
tionality between shear rate and segregation rate, and 
sets the overall duration of the process. This segrega- 
tion rate has previously been observed to vary with such 
parameters as particle size-ratio and confining pressure 
[15] , We find that the model successfully captures the 
existence of a fast mixing process followed by a slower re- 
segregation process. However, the model re-segregates in 
finite time and does not replicate the observed behavior 
from the experiment, where a completely re-segregated 
state is approached only exponentially in time. 
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FIG. 2: (a) H(t) for 11 experimental runs (gray lines), with 
thicker line representing the run for which we took movies 
and measured the velocity profile. The thick black line is the 
average of the 11 runs, used for comparison with the model 
in Fig. [7] Positions of particles (o indicates small, x indicates 
large) overlaid on images taken through the window in the 
outer wall (b) initially, t = sec, (c) in a mixed state, t = 
127 sec, and (d) in the re-segregated state, t = 1820 sec. 



II. EXPERIMENT 

The experimental apparatus is an annular Coucttc cell 
filled with a bidisperse mixture of spherical glass particles 
confined by cylindrical walls at inner radius (25.5 ± 0.1) 
cm and outer radius (29.3 ±0.1) cm. We apply shear via 
a circular bottom plate rotating at a constant frequency 
/ = 49 ± 0.5 mHz, approximately 3 rpm. A heavy top 
plate sits within the annular gap and is free to move 
vertically, but is partially suspended by springs to reduce 
the pressure it applies. The compressive force P applied 
to the particles is (0.36±0.008) mg, where mg is the total 
weight of the particles and the variation in force is due 
to the stretching of the springs. The top and bottom 
plates have a rubberized surface to increase friction with 
the particles, while the stationary cylindrical side walls 
are constructed of aluminum. 

At the start of each experiment, we prepare a flat layer 
of 2.0 kg of large particles (diameter di = 6 mm) at the 
bottom of the annulus, followed by a flat layer of 2.0 kg 
of small particles (ds = 3 mm). This configuration is 
shown schematically in Fig. [l] We consolidate the layers 
by compression prior to the beginning of the run, with 
the average initial height H(0) = 41.2 mm. 

Fig. [2] shows the height H(t) of the top cell as a func- 
tion of time for several different runs under the nominally 
identical conditions. While the variability from run to 
run is considerable (and typical of granular materials), 
the important features are common to all of the runs. 
From the initial normally-graded configuration, the ap- 
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plied shear causes the lowest (large) particles to move 
horizontally due to friction with the rubberized bottom 
plate. As shear begins, the material must dilate in order 
to deform (Reynolds dilatancy). After this initial rapid 
expansion, H(t) provides information about the progres- 
sion of mixing/segregation in the aggregate. Because the 
mixing process causes small particles to partially fill the 
voids between the large particles, the total height of the 
aggregate decreases. For all runs, the height reaches a 
minimum H m i n at t w 80 sec, with AH ss 1 mm. Fur- 
ther shearing serves to re-segregate the particles, with the 
large particles ending at the top of the cell and the small 
particles at the bottom. During this process, H rises 
back to a height greater than its initial value. This re- 
segregation process would in theory continue until there 
were no longer a mixture of particles in the central re- 
gion. In practice, however, the particles do not fully 
re-segregate and a few large particles remain within the 
lower layer, as can be seen in Fig. [2ji. 

To measure the velocity of the particles, we observe 
the outer layer through a window of approximate width 
10 cm, using a digital video camera operating at 450 
Hz. We recorded images during three time intervals, 
each with an approximate duration of 10 minutes (a to- 
tal of around 10 6 images), separated by intervals of sim- 
ilar duration during which images were transferred from 
the camera to the computer. The system reaches a re- 
segregated state after approximately tf = 700 seconds. 

From each image, we first identify the center of each 
particle by convolution with a circular kernel chosen to 
match either ds or cLl- We perform the convolution 
twice (once for each particle-diameter) and then screen 
for mis-detections and double-detections. Fig. [2] shows 
particle centers for (b) particles in the initial, normally- 
graded configuration, (c) the mixed state, and (d) the 
re-segregated, inversely-graded state. Because the par- 
ticular configuration visualized at the wall only measures 
the state of a small portion of the system, we use the im- 
ages only for calculating the velocity profile, and use H(t) 
to probe the average degree of mixing/segregation. 

To characterize the shear, we are primarily interested 
in the average horizontal velocity u and how it depends 
on depth, time, and particle size. We assemble parti- 
cle trajectories from the list of particle positions associ- 
ated with each image, considering each of the two sizes 
separately. For each trajectory, we calculate the instan- 
taneous velocity of the particle from the slope of a lin- 
ear fit over a duration appropriate to the average speed 
of the layer. This analysis was repeated separately for 
each of the three 10 minute intervals to check for time- 
independence; we include all images with t < 3100 sec- 
onds in order to improve our statistics. We observed a 
steady-state velocity profile for all three intervals (shown 
in Fig. [3]), after an initial sa 37 sec transient which is ex- 
cluded from the analysis. We observered that u(z) was 
approximately the same for both large and small parti- 
cles. 

We measure the average horizontal velocity as a func- 
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FIG. 3: Measured velocity profile u{Zi) (•) for (a) full cell 
height, with boundary layers of thickness above and below 
the dashed horizontal lines and (b) within the central region, 
scaled so that z £ [0,1]. The velocities are scaled by u = 
5.5 mm/sec so that u(0) = 1. The solid curve is the fit to 
Eq. |TJ. Bars represent width of the velocity distribution at 
half the height of the peak. 



tion of depth by dividing the ensemble of trajectories into 
a discrete set of bins centered at positions Zi- Within 
each bin, we plot a probability distribution of the veloci- 
ties and fit a parabola to the peak. This maximally-likely 
velocity, u(zi), is plotted in Fig. [3] with the width of the 
distribution at half the height of the peak represented by 
the horizontal bars for each Zi. Note that the velocity 
fluctuations ("granular temperature") are largest within 
the shear band, and that thermal gradients have been as- 
sociated with segregation in kinetic theory models [HI [TO] ■ 
Since the segregation is driven by shear in our model, 
we calculate the experimental shear rate 7(2) = \du/dz\ 
from velocity data by finite-differences, shown in Fig. [4] 
Note that the bottom and top layers, of depth approxi- 
mately dz, , show dramatically higher shear rates than the 
bulk. In the figures, we mark the boundaries between the 
bulk and the two layers by horizontal dashed lines. In or- 
der to model the long-term mixing and re-segregation, we 
focus on the central region of the flow and scale the z axis 
so that z = at the bottom of this region and z = 1 at 
the top. Within this region, we observe that the velocity 
profile is well-described by an exponential of the form 

u(z) = it e~ 2/A + c, (1) 



4 




FIG. 4: Dimensionless shear rate 7 = \du/dz\ (•) for (a) full 
cell height, with boundary layers above and below the dashed 
horizontal lines and (b) within the region 2 G [0, 1]. Solid line 
is the fit from Fig. [3) plotted as Eq. (2|. 



where A is related to the width of the shear band, uo + c 
is the velocity at z = 0, and c is a constant representing 
solid-body rotation. The fit to this form is plotted as a 
solid line in Fig. JSJd, scaled so that u(0) = 1. The fit 
parameters are A = 0.205, uo = 0.866, and c = 0.134. 

To obtain 7(2) for use in the model, we differentiate 
u(z) in Eq. ([IJ of u(z) and find 



7(2) 



Uo 

A 



-z/\ 



This function is plotted in Fig. |4Jd for comparison with 
the experimentally-determined profile obtained by finite- 
differencing. In we will relate the local segregation 
rate to the parameters uq and A. 



III. SEGREGATION MODEL 

To simulate the segregation and mixing in the exper- 
iment, we begin with an equation for the conservation 
of mass and impose our modeling assumptions to mimic 
conditions in the experiment. In the model, we consider 
the kinetic sieving mechanism but not diffusion. The 
conservation of mass equation is |13j : 
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dt dx 



d_ 

dz 



d 



+ — ((j)u) + — (cj)v) + — (^w) + — (q{z)cj)(cj)-l)) 



= 0, 
(3) 



where = 4>(z, t) is the concentration of small particles 
at position z and time t. The velocity (it, v, w) is the bulk 
velocity of small particles. Kinetic sieving is modeled by 
a modification w s — q{4>~ 1) to the vertical component w 
of bulk velocity, in which q = q(z) is the local segregation 
rate. 

We assume that the components ii,kj of the bulk ve- 
locity (in the y, z directions, respectively) are negligible. 
That is, we assume there is essentially no motion across 
the flow, and that the vertical component of velocity of 
small particles is dominated by the effect of segregation. 
In the annular Couette geometry of the experiment, the 
flow is uniform in the ^-direction (the angular direction) ; 
we assume u = u(z) is independent of horizontal position 
(x, y) and time t. Finally, we assume that the segregation 
rate q{z) depends only on the vertical variable z and is 
proportional to the shear rate j(z): 



q(z) = sj(z) = q e z/X , 



9o 



u 



(4) 



The constant of proportionality s sets the time scale in 
the model and represents the rate at which mixing and 
segregation processes occur. In previous experiments, we 
observed that the timescale of both the mixing and seg- 
regation processes is a function of particle size ratio and 
confining pressure |18j . 

With these assumptions neglecting small velocities, the 
model ([3]) reduces to the scalar conservation law 







+ -(g(z) M-1))=0. 



dt d 



(5) 



We set an initial condition corresponding to the begin- 
ning experimental configuration of a layer of large parti- 
cles above a layer of small particles: 



0(*,O) 




(2) with zq= h. Boundary conditions 



#0,t) = 1, 0(1,*) =0 



(6) 



(7) 



ensure that there is no flux of the particles through the 
upper and lower boundaries. 

In May et al. [14] , we constructed the solution of Eq. ^ 
[7]) using characteristics and shock waves. The early time 
solution is a rarefaction wave centered at z = zq, t — 
0, that represents the initial mixing, in which <j) varies 
continuously from <j) = to <f> = 1. The solution in the 
rarefaction fan is characterized by a pair of equations 



<f>(z,t) = -<r o *0 o (l - o )e-" o/A + O 



(8a) 



e z/X - e 2 °/ A = -0 O (1 - o )e- 2 °/ A (a o t) 2 + (20 o - l)a t, 

(8b) 

where a = q / X. In these equations, < 4>o < 1 labels 
a specific characteristic at * = 0, z = zq; cf>o(z,t) can be 
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FIG. 5: Numerical solution of the initial value problem Eq. |5] 
-[7]), with parameters Uo = 0.866, A = 0.205, zq = §. For this 
plot, we set the segregation rate s = 16.2 s _1 so that tf — 1 
(non-dimensional) at the final, fully-segregated state. 



In this section, we relate the mixing and re-segregation 
stages in the height evolution to predictions of the model 
([5]), shown in Fig. [5} 

The solution cf>(z,t) of the model is specified at each 
point z in the fixed domain < z < 1, and at each time 
t. However, the physical domain varies with time, due 
to the changes in height H(t). To relate the evolution of 
cf>(z, t) to H(t), we calculate a proxy height that takes into 
account the different packing densities of the material at 
different concentrations <j) of small particles. The local 
packing density p is the ratio of the solid local volume 
occupied by the mixture of particles to the total local 
volume, i.e., including void space. The total volume of 
the annular region is the product of the cross sectional 
area A and the position H of the top plate. For a local 
solid volume of particles dV p across a cross section of the 
annulus, in a horizontal layer of small (physical) thickness 
dh, we write 



found by solving the quadratic equation ( 8b ) , and choos- 
ing the relevant solution. Then Eq. (8a I is an explicit 
formula for cj)(z,t). 

The rarefaction reaches the bottom and top boundaries 
when the characteristics (j> = 1,0 = reach z = 0, z = 1, 
respectively. This corresponds to the first small parti- 
cle reaching the bottom plate and the first large particle 
reaching the top plate. These events occur at times 



to = -(e Za/X - 1), 
9o 



h = -(e 1/A -e z ° /A ), (9) 



respectively. Subsequently a layer of small particles 
grows from z = 0, and a layer of large particles grows 
from z — 1. The interfaces z — T (t),z = Ti(t) between 
the rarefaction and these layers are shock wave solutions 
of the conservation law ([5]), and consequently evolve ac- 
cording to the Rankinc-Hugoniot condition, which for 
Eq. (5} is a differential equation for each shock: 

T' (t) = g(r o (i))0(r o (t), t), t > t , r (t ) = 0, (10a) 

ri(t) = ^r^x^r^),*) - 1), t> t x , T 1 {t 1 ) = 1. 

(10b) 

When these shocks meet, re-segregation is complete, and 
the solution consists of a stationary shock separating the 
upper layer of large particles from the lower layer of small 
particles. By mass conservation, the position of the in- 
terface is z = 1 — zq. A contour plot of the solution 
corresponding to parameter values calculated from the 
experimental data is shown in Fig. [5j 



IV. COMPARISON OF MODEL AND 
EXPERIMENT 

In the experiment, the height H{i) initially increases 
due to Reynolds dilatancy, then decreases sharply as the 
material mixes (see Fig. |2^). Subsequently, H(t) in- 
creases at a slower rate as the material re-segregates. 



dV p 
Adh' 



(11) 



We will assume that the local packing density p depends 
only on the local concentration <fi(z,t). That is, we take 
p = p{4>), independent of other factors such as the shear 
rate. Thus, for each dz (the height element in the model 
domain, < z < 1), there will be a local contribution 
to H{t) which depends on the local <f>. In order to con- 
vert from the model domain to the physical domain, we 
integrate to obtain a proxy height H{t): 



H(t) 



H(t) 



dh 



A 



dz 



(12) 



where V p is the total solid volume of the granular mate- 



rial. We note that in Eq. (12 1, -f has units of length, 



and the integral over dz is nondimensional. Another way 
to interpret this equation is that the integral calculates 
the height of the sample in the mathematical model and 
the coefficient -£■ relates the model height H(t) to the 
physical height H(t). 

To determine an appropriate function p((f>), we first 
consider the monodisperse cases in which <j> = or <f> = 1 . 
In these cases, the lower limit for p is known as random 
loose packing, (p RLP w 0.55) and the upper limit as ran- 
dom close packing (p RCP w 0.64) [Ti?l 120) . However, for 
bidisperse mixtures, the packing density depends on the 
relative composition of the mixture. Since small parti- 
cles can fit within the spaces between large particles, the 
packing density is larger for a bidisperse mixture than 
for a monodisperse sample. Indeed, the maximum pack- 
ing density for a bidisperse mixture with this size ratio 
has been observed to be around p RCP > bl = 0.67 to 0.69, 
depending on the method used [TT1 |2"T1 122"] . 

Data from several numerical and experimental stud- 
ies [17l [23j [24] of the random close packing of a static 
bidisperse mixture of spheres supports an approximately 
triangular shape for the concentration map p(<fi), with a 
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FIG. 6: Contour plot of H m i n — H m in (solid lines) and AH — 
AH (dashed). Symbol • marks the intersection of the two zero 
contours at V p / (Ap m i n ) = 42.0 mm and Ap/p m i n = 0.072. 



FIG. 7: Comparison of the experimentally-measured H(t) 
(solid line, red, see Fig. [5| and the calculated proxy height 
H (dashed line, black) for the parameters given in Fig. Im 



maximum at <b r — 0.275. The function 



p(<j>(z,t)) 



< 



> 4>c 



(13) 



falls off, by an amount Ap, to minimum packing density 
Pmin at a 4> = 0, 1. Kristiansen et al. |17j report a value 
of Pmin = 0.628, which is below both the monodisperse 
and bidisperse p RCP . 

Note, however, that there is a significant difference be- 
tween the conditions under which these studies deter- 
mined the packing density, and the conditions of our ex- 
periments: our measurements take place in a sheared 
system rather than a static one. The packing density for 
sheared granular materials is typically less than for static 
packings, due to Reynolds dilatency. However, we are not 
aware of any measurements, analogous to those reported 
in Kristiansen et al. |17| , for p((f>) in sheared granular ma- 
terials. (One can imagine that segregation in fact makes 
such measurements rather challenging.) Therefore, we 
must estimate the values of p rn in an d Ap in the concen- 
tration map (13 1 for our bidisperse mixture of granular 



materials under shear. We expect that the value of p m in 
may be less than the value found in Kristiansen et al. |17j 
due to the shear, but is probably not less than the ran- 
dom loose packing limit for monodisperse spheres, p RLP . 
Therefore, we consider 0.55 < p m i n < 0.64 as a reason- 
able range of values. Typical reported values of Ap range 
from 0.025 to 0.063 [17] ■ Since p RCP > bl provides a larger 
upper limit for bidisperse mixtures, we consider a large 
range of values, 0.02 < Ap < 0.14. 

Substituting Eq. (13 1 into (12), we observe that there 
are only two free parameters: A J P , sets the overall 
height of the system and sets the amount of com- 
paction/expansion. From measurements of apparatus di- 
mensions and particle sizes, we find 5f = 24.2 ±1.2 mm. 
This range, together with the ranges for p and Ap 
given above, predict that we should consider parameters 



the values which best capture the compaction and ex- 
pansion process, we perform the integral ( 12 1 for pairs of 

( 



Ap 



) values, and determine the resulting H„ 



37 < 



Ap„ 



< 44 and 0.07 < 



< 0.26. To select 



Ap min 'j>„..,. 

and AH. In Fig. |6| we show contour plots of the 
difference between the experimentally-measured values 
and the proxy-calculated value over this full parameters 
range. The best parameter choice lies at the intersection 
of the two zero-contour lines. We find these values to be 
A -^— = 42.0 and = 0.072. The first parameter is 
to be expected given the height of the system in Fig. [2] 
and the second is at the lower end of the expected range 
of values. 

A third free parameter is the segregation rate s in 
Eq. (|5| that sets the time scale for the dynamics. From 
the solution of the initial value problem, we find that 
the choice s = 0.023 predicts the ending time for the 
solution equal to the experimentally-determined time 
tf = 700 sec. Fig. [JJ provides a direct comparison of the 
experiment and the model over this time interval. We ob- 
serve that the proxy height decreases more rapidly than 
it climbs, in agreement with the experimentally observed 
H(t); this feature arises without including separate mix- 
ing and segregation rates in the model. The microme- 
chanical origins of this property remain to be investigated 
in future work. In addition, the rate at which it decreases 
is close to the observed rate. At later times, the shapes of 
H{t) and H (t) no longer agree, with H(t) exhibiting both 
a flatter minimum and a faster re-segregation rate. It is 
particularly notable that the model solution predicts re- 
segregation in finite time, in contrast to the exponential 
approach to re-segregation observed in the experiment 



V. CONCLUSIONS 

We have performed quantitative experiments to com- 
pare predictions for a recent extension |14j of the Gray- 
Thornton model of granular segregation |13j to a flow in 
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which the shear rate is non-constant. This comparison 
additionally requires the construction of a model connect- 
ing the local concentration of small and large particles 
to changes in the local packing density (the concentra- 
tion map). When these two elements are combined, we 
can compare the temporal evolution of the height of the 
segregating system to a proxy height calculated from the 
continuum model. Several features of the experimentally- 
measured height dynamics H (i) are well-captured by the 
proxy height: (1) a fast mixing timescale is followed by a 
much slower re-segregation timescale, (2) we can model 
the compaction and re-expansion process using reason- 
able parameter values in the concentration map, and (3) 
the slope of the height curve during the mixing phase is 
in approximate agreement with that of the model. 

Some prominent features are missing from the model: 
there is no means to account for Reynolds dilatency, 
and the model segregates in finite time rather than 
exponentially-approaching a final, re-segregated state. 
This latter point of disagreement shows the limitation of 
using a continuum mixture-theory model for a discrete 
process, especially where the number of discrete objects 
(the particles) is comparatively small. In the experiment, 
it is easy to see how this finite-size effect takes over af- 
ter the continuum model predicts complete segregation. 
Once most of the large particles have reached the upper 
layer, it becomes increasingly difficult for the few large 
particles remaining in the lower region containing mostly 
small particles to segregate to the upper layer. When 



large particles do not have large particle neighbors, the 
small particles cannot sieve between the large particles 
and thus move the large particles upward. In fact, some 
of the large particles never make it to the upper layer, 
even after runs of several days duration, instead remain- 
ing trapped among the small particles. 

Despite the differences between the experimental and 
proxy heights, the proxy height calculated from this sim- 
ple continuum model captures the qualitative features of 
the experimental height time series. We have extended 
the Gray-Thornton model, which was developed for the 
case of an avalanche (uniform shear), to model nonuni- 
form shear, as occurs in a one dimensional configuration 
in an annular Couette cell. It is noteworthy that a con- 
tinuum model applied to a small scale granular system 
successfully captures the main phenomena of mixing and 
re-segregation. 
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